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ABSTRACT 

We present the i?-band luminosity function for a sample of 18678 galaxies, 
with average redshift z — 0.1, from the Las Campanas Redshift Survey (LCRS). 
The luminosity function may be fit by a Schechter function with M* = —20.29 ± 
0.02 + 51og/i, a = -0.70 ± 0.05, and 0* = 0.019 ± 0.001 Mpc-^ for absolute 
magnitudes -23.0 < M - 51og/i < -17.5 and h = i/o/(100 km s"^ Mpc'^). 
Over the same absolute magnitude range, the mean galaxy density is 0.029 ± 
0.002 Mpc~^ for a volume extending to cz = 60000 km s~^. 

We compare our luminosity function to that from other redshift surveys; in 
particular our luminosity function normalization is consistent with that of the 
Stromlo-APM survey, and is therefore a factor of two below the normahzation 
implied by the bj ~ 20 bright galaxy counts. Our normalization thus indicates 
that much more evolution is needed to match the faint galaxy count data, com- 
pared to minimal evolution models which normalize at bj ^ 20. Also, we show 
that our faint-end slope a = —0.7, though "shallower" than typical previous 
values a = —1, results primarily from fitting the detailed shape of the LCRS 
luminosity function, rather than from any absence of intrinsically faint galaxies 
from our survey. 

Finally, we find that the faint end of the luminosity function is dominated by 
galaxies with emission lines. Using [Oil] 3727 equivalent width W\ = 5 A as the 
dividing line, we find significant differences in the luminosity functions of emis- 
sion and non-emission galaxies, particularly in their a values; emission galaxies 
have Schechter parameters M* = —20.03 ± 0.03 + 51og/i and a = —0.9 ± 0.1, 
while non-emission galaxies are described by M* = —20.22 ± 0.02 + 51og/i and 
« = — 0.3±0.1. The average [Oil] 3727 equivalent widths do not change signifi- 
cantly with redshift, consistent with a star formation rate which stays constant 
over the depths sampled by the LCRS. This result holds for galaxies of different 
luminosities, and over the respective redshift ranges that these galaxies may be 
observed, in particular up to about z — 0.2 for galaxies brighter than M*. 

Subject headings: cosmology: observations — galaxies: distances and redshifts 
— galaxies: luminosity function, mass function — surveys 
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1. Introduction 

The galaxy luminosity function is one of the most basic and fundamentally important 
quantities in observational cosmology. For example, it is essential for the analysis of large- 
scale structure statistics from redshift surveys, and knowledge of the local luminosity function 
is necessary for the proper interpretation of the faint galaxy count data and galaxy evolution 



at high redshift (e.g., |Koo fc Kron 19921) 



Though much work has gone into determination of the local z = optical luminosity 
function (see [I'elten 1977| ; [Kirshner et al. 19"79| , 1983; [Efstathiou et al. 1988| ; [Binggeli et 



|al. 1988| ; and recent work by [Loveday et al. 1992| ; [Marzke et al. 1994a| ), there is still some 



uncertainty in its faint end as well as in its normalization (e.g., |Colless 19951 ). In this paper 
we investigate the local luminosity function by presenting results for a sample of 18678 Las 
Campanas Redshift Survey (LCRS) galaxies, selected from CCD photometry in a "hybrid" 
red Kron-Cousins i?-band (see § ^ below). This is the largest galaxy sample for which the 
luminosity function and mean density have been computed. The combination of the mean 
survey depth of 2; = 0.1 and the large survey volume (1.0 x 10^ Mpc^ out to cz = 60000 
km s~^) imply that our results are not biased by local inhomogeneities. In addition, the 
spectra of LCRS galaxies allow us to study the luminosity function for samples divided by 
[Oil] 3727 emission strength, a rough indicator of the star-formation rate ( [Kennicutt 1992| ). 

We describe relevant aspects of the LCRS construction in § We then detail our 
methods for measuring the galaxy luminosity function and mean density in § 0. Our results 
are presented in § ^, and we compare them to those of other redshift surveys in § ^. Some 
implications of our results for faint galaxy counts and galaxy evolution are discussed in § ^ 
We summarize our conclusions in § ^ 



2. The Las Campanas Survey Data 

We begin by describing the Las Campanas survey parameters; detailed descriptions 
of the LCRS may be found in Shectman et al. (1992, 1995, 1996), Tucker (1994), and 
Lin (1995). There are 23690 galaxies, with an average redshift z = 0.1, which lie within 
the survey photometric limits and geometric borders considered in this paper. The survey 
galaxies were selected from a CCD photometric catalog obtained from drift scans taken 
on the Las Campanas 1.0 m Swope telescope with a Gunn r filter (P?huan fc Gunn 1976| ). 
Follow-up spectroscopy was obtained using the multi-object fiber-optic spectrograph on the 
Las Campanas 2.5 m DuPont telescope. The survey geometry is that of six "slices," 1.5° in 
declination by 80° in right ascension, three each in the North and South galactic caps. 
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The survey data may be divided into two main parts: the first 20% of tlie data was 
obtained using a 50-object fiber system, while the remaining 80% of the data was taken with 
a 112-object system. The survey shces were built up by observing 1.5° by 1.5° spectroscopic 
fields, with a maximum of 50 or 112 galaxies observed per field. The survey photometric 
limits were chosen so that there would be typically more targets per field than available fibers, 
and we selected targets at random among those that met the selection criteria. Because we 
generally do not re-observe any of our fields, we keep track of the variable field-to-field 
sampling fractions fi in our subsequent statistical analyses. 

To be precise, though our photometry was obtained through a Gunn r filter, the calibra- 
tion was done relative to standard stars ( [Graham 1981| ; [Graham 1982| ) in the Kron-Cousins 



i?-band. Our magnitudes are thus in a "hybrid" Kron-Cousins R system. Nevertheless, the 
difference between our hybrid Kron-Cousins R system and the true one is small; compari- 
son of photometry of LCRS galaxies obtained in both systems yields a zero-point difference 
of < 0.1 mag and a color term of < 0.2 mag per magnitude of the color (Johnson V— 
Kron-Cousins R). See Tucker (1994) and Shectman et al. (1996) for more details. 

Objects were selected for the survey using isophotal magnitudes m, applying both faint 
and bright isophotal magnitude limits. The limiting isophote is 15% of sky, or approximately 
23 mag arcsec"^. In addition, a central magnitude rric, approximately measuring flux that 
would enter through the 3.5"-diameter aperture of a fiber, is also computed for each object. 
An additional cut line in the m-rric plane is then imposed in order to exclude low central 
surface brightness objects, for which redshifts are more difficult to obtain. For the 50-fiber 
data, the nominal isophotal magnitude limits are 16.0 < m < 17.3, and the central surface 
brightness cut excludes about 20% of objects otherwise meeting the isophotal limits. For 
the 112-fiber data, the nominal isophotal limits are 15.0 < m < 17.7, with exclusion of just 
the lowest 4-9% of objects by surface brightness. We illustrate our isophotal and central 
magnitude selection criteria for a sample of LCRS galaxies in Figure |I]. There are small 
field-to-field variations from the nominal photometric limits; these are accounted for in our 
analyses. We will find it convenient to divide our data into four subsets by galactic cap. 
North and South, and by fiber number, 50 and 112 fibers; we will refer to these subsets as 
N50, S50, N112, and S112. Table ^ summarizes details concerning these four samples. Also, 
it turn outs that our fiducial sample will be the combined North plus South 112-fiber data 
set, which we will refer to as NS112. 

Throughout this paper we adopt a Hubble constant Hq = 100 h km s^^ Mpc^^ and a 
deceleration parameter qq = 0.5. Where the h dependence is not explicitly indicated, we 
have used h = 1. We also assume a uniform Hubble flow so that the absolute magnitude M 
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of a galaxy of isophotal magnitude m and redshift z is given by 



M = m 

where the luminosity distance di is 

c 



25 -5 log dL - k{z) 



di 



qoz + (1 - go)(l - \/l + 



and where the /c-correction k{z) is taken to be 



k{z) = 2.51og(l + 2) 



(1) 



(2) 



(3) 



as appropriate in the mean for an r-band-selected galaxy sample ( [Tucker 1994 ; [Kirshner 



et al. 19831) . Note that galaxy pecuhar velocities, of order ai2/V2 ~ 350 km s ^ (cri2 is 



the pairwise peculiar velocity dispersion; see [Lin et al. 1996E| ; [Marzke et al. 1995[ ), and 
our velocity measurement errors, on average 67 km s~^ ( P^in 1995 ), are both quite small 
compared to the mean survey redshift cz = 30000 km s~^. We therefore neglect the effect of 
peculiar velocities and velocity errors in the calculation of absolute magnitudes. Similarly, 
the effect of using a different value of go is also small: M(go') — M{qo = 0.5) ~ (go' — 0.5)z, 
e.g., —0.04 magnitude at our average z = 0.1 for the case go' = 0.1. 



3. Methods 

3.1. Luminosity Function 

We compute the luminosity function using two related methods which are unbiased by 
density inhomogeneities in the galaxy distribution. These are the parametric maximum- 
likelihood method of Sandage, Tammann, & Yahil (1979; hereafter STY) and the non- 
parametric stepwise maximum-likelihood (SWML) method of Efstathiou, Ellis, & Peterson 
(1988; hereafter EEP). We describe these methods below, and we will discuss some modifi- 



cations needed in the analysis of the LCRS sample in § p]2 



Consider a galaxy i observed at redshift Zi in a flux-limited redshift survey like the 
LCRS. Let mmin,i and mmax.i denote the apparent magnitude limits of the field in which 
galaxy i is located. Let (/>(M) be the differential galaxy luminosity function which we want 
to determine. Then the probability that galaxy i has absolute magnitude Mj is given by 

/ 0(M)dM , (4) 
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where the absolute magnitude hmits at Zi are 

1 -^^max ( ■^^i ) J 1 '^maxji 



2^-^\ogdLi-k{zi) . (5) 



The luminosity distance di and the k correction k{z) for LCRS galaxies are computed as 
described in § ^. We can then form a likelihood function £ for having a survey of galaxies, 
with respective absolute magnitudes Mj, by multiplying the probabilities pi. 

N 

£ = p{Mi,...,Mn\zi,...,zn) = Y[pi , (6) 

i=l 

or 

ln£ = ^nn0(Mi) -In / (p{M)dM\. (7) 

In the STY method one assumes a parametric model for 0(M), and the parameters describing 
0(M) are determined by maximizing the likelihood £, or equivalently ln£, with respect 
to those parameters. In our case we take as our model for 0(M) the Schechter function 
( [Schechter 197^ ) 

0(M) = (0.4 In 10) 0* [ioO-4(A^*-Af)]i+a exp[-10°-"(^*-^^)] , (8) 

and use the STY method to find the characteristic magnitude M* and the faint-end slope a. 
The normalization 0* drops out in equation (^) and has to be determined using a different 
method to be described below in § p.3| . Also, error ellipses in the M*-a plane may be drawn 
by finding the contour corresponding to 

ln£ = ln£max-|Ax^ (9) 

where Ax^ is the change in appropriate for the desired confidence level and a distri- 
bution with 2 degrees of freedom. 

Alternatively, one does not have to assume a particular functional form for 0(M). 
Rather, in the EEP stepwise maximum likelihood method (SWML), the luminosity func- 
tion is taken to be a series of Np steps, each of width AM in absolute magnitude: 

0(M) = 0fc , Mfc - AM/2 <M <Mk + AM/2 , k = I, . . . , Np . (10) 

In this case equation (0) may be rewritten as 

N Np 

ln£ = ^^1^(M,-Mfc) ln0fc 

i=l k=l 

- Y.^n\j2(PkAM H[Mk,M^,^{z,),M^U^,)] \ , (11) 

1=1 I k=l I 



- 7- 



where the functions W and H are given by 

Tx./,. f 1, Mk- AM/2 < Mi < Mk + AM/2 

otherwise ' ^^^^ 

i^■[Mfc,M^i„(^^),M^ax(^^)] = 

min[Mfc + AM/2, M^^axl^i)] - max[Mfc - AM/2, M^i^{zi)] , 

if Mfc + AM/2 > Mminizi) and M^ - AM/2 < Mmax(-2i) • (13) 
0, otherwise 



We then maximize ln£, equation ([TT|), with respect to the parameters 0^, and subsequently 
solve for them iteratively using 

4>kAM = |^iy(M,-M,)| ^ (14) 

I H[Mk, M^ir^iZi), M^ax(^i)] / E H[Mj, MrainiZi), M^,,{Zi)] 

i=l [ / j=l 

As in the STY method, the normalization of the 0^ is not determined, so we introduce a 
constraint similar to that given by EEP to fix the normalization: 

g = Y.(l)kAM V{Mk)/V{Mj) -1 = 0, (15) 

k=l 

where V{M) gives the volume (including relativistic effects) over which a galaxy of absolute 
magnitude M can be seen in a flux-limited survey, and Mj is a flducial absolute magnitude, 
which we take to be —20 ~ M* as appropriate for the LCRS. Also following EEP, we 
can estimate the variances of the parameters (j)k from the first Np diagonal elements of the 
covariance matrix given by 

coY{<Pk) = r\<f)k) , (16) 

where l{4>k) is the information matrix 



92 In £/d(j)id(j)j + {dg/d(j)i){dg/d(j)j) dg/d(j)j 
dg/d(l>i 



(17) 

k,4>j=4'k 



We can test the goodness of fit of the STY solution by comparing it to the SWML 
solution using a likelihood ratio test described by EEP. Let £2 be the likelihood for the SWML 
solution given by equation (pT]), and let £1 be the likelihood, again from equation (pUD, but 
with 0fc set to 0sTY,fc, where 



^STY,fc 



(Pl^^{M)V{M)dM / / 0sTY(M)F(M)ciM 

A4-AA//2 / JMk-AM/2 
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As discussed in EEP, here we use this definition of £i for the STY solution, rather than equa- 
tion in order to match the binning present in equation (|11]) but absent from equation (J^). 
EEP show that the hkehhood ratio 21n(£i/£2) approximately follows a distribution with 
Np ~ 3 degrees of freedom, giving us a means for testing the quality of the STY solution. 

Alternatively, we can also check the goodness of the STY fit by comparing the observed 
absolute magnitude distribution dN{M) against the distribution dNsTY^M) expected from 
the STY solution and the observed redshift distribution of the survey galaxies (|Yahil et] 
al. 19911) : 

dNsTviM) = Y,P{M\zi)dM = ^ <^ (p{M)dM / / (p{M)dM \ . (19) 

i = l i = l I / •'^Iminiz^) J 

A test can then be used to assess how well the STY prediction matches the observations. 



3.2. Modifications 

For thoroughness in the analysis of the LCRS data, we need to consider four additional 
details which will modify the basic methods described above to calculate the luminosity 
function. The most important effect is that of central surface brightness selection, while the 
others make fairly small changes in 0(M). 

(1) Magnitude Errors. From repeat observations we find that our rms isophotal mag- 
nitude error is a = 0.1 magnitude. We make the rough approximation that the magnitude 
errors are Gaussian distributed with dispersion a, so that the observed luminosity function 
(pobs is a convolution of the true luminosity function with a Gaussian magnitude error 
distribution (EEP). Thus 

1 



^obs 



27TCX 



c^{M')e-^^''-^'^ dM' , (20) 



and we use (pots in place of in the parametric STY equations. For simplicity, we will 
not make this substitution in the non-parametric SWML equations. Accounting for this 
dispersion of 0.1 magnitude makes small systematic changes in the best-fit values of the 
Schechter parameters for our sample: AM* ^ +0.03 and Aa ^ +0.03. 

(2) Variable Field Sampling Fractions. To account for this effect we weight galaxy z's 
entry in the likelihood equations and (|ll|) by the inverse of the sampling fraction fi of 
the field containing that galaxy ( |Zucca et al. 1994]) : 

Inp, . (21) 

ft 
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Including this effect changes the Schechter parameters very httle: both |AM*| and |Aq;| are 
< 0.01. 

(3) Apparent Magnitude and Surface Brightness Incompletenesses. Our success rate in 
identifying spectra as galaxies and stars drops slightly at fainter apparent magnitudes and 
also drops slightly close to the limiting central surface brightness cut line (Figure |l|). We 
define a completeness factor F(m, mc — rnc^umit)-, which is a function of isophotal magnitude 
m and of "distance" — rni.iimit{m) from the central magnitude cut line. We take F to 
be the fraction of spectra we identified as galaxies or stars in each bin on a grid of m and 
rric — mc^iimitijn) values. In Figure 0, we plot the two orthogonal projections Fi{m) and 
^2(^0 — fnc,iimit) of the full completeness function F for each of the N50, S50, N112, and 
S112 samples. From the figure we note that the completeness does not vary much with either 
apparent magnitude or surface brightness; the completeness only drops to about 0.95 at the 
faint isophotal magnitude limits or at the central surface brightness limits. To account for 
this effect, we weight galaxy z's entry in the STY and SWML likelihood equations by the 
inverse of the value of F appropriate for that galaxy's isophotal and central magnitudes: 

Inp, — ^ . (22) 

F{mi,mci - rriciimiti) 

The change introduced is negligible for M*, |AM*| < 0.01, but a does change systematically 
by a small amount, Aa ~ —0.04. (The greater incompleteness at fainter m translates to a 
bias against intrinsically fainter galaxies and hence a more positive a, if left uncorrected.) 

(4) Central Surface Brightness Selection. We empirically find in the LCRS a correlation 
between a galaxy's absolute luminosity and that galaxy's distance from our central surface 
brightness cut line. We demonstrate this by culling out galaxies from our N112 subset to 
match the more restrictive surface brightness cuts of the N50/S50 and S112 subsets (§ ^ 
and Table p. Figure |^ shows that the effect of the more stringent surface brightness cuts is 
to reduce the sampling completeness G{M) as the absolute magnitude M becomes fainter, 
where G{M) is simply the fraction of galaxies at M which remain after the low-surface- 
brightness galaxies have been removed. For the N50, S50, and SI 12 samples, we can account 
for this bias against faint galaxies by weighting galaxy z's entry in the likelihood equations 
by the inverse of G at that galaxy's absolute magnitude Mf. 

We thus correct the surface brightness sampling of the N50, S50, and S112 subsets to that 
of the N112 subset. For the 50-fiber data samples, AM* ^ —0.04 and Aa ^ —0.15; for the 
S112 sample, AM* = —0.01 and Aa = —0.06. We do not attempt to further extrapolate 
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and correct for the N112 sample's own central surface brightness cut, though we expect the 
change to be |Aa| < 0.1 from the example of the SI 12 data set. 

We take account of the latter three effects (2-4) together by the replacement 

Inp, — . lnp„ = -— VrTTT^ ' (^4) 

that is, we weight by the inverse of the product of the three completeness factors. However, 
because the weighting factors Wi are greater than one, the likelihood ln£ will be larger in 
absolute value than when no weightings are applied. The likelihood will thus resemble that 
of a larger sample, leading to systematic underestimates in the confidence ellipses and errors 
we compute (|Glazebrook et al. 1995c|) . To remedy this, we renormalize the weights Wi so 



that Y^f^i Wi = N, where N is the number of galaxies in the sample of interest. The original 
weights are multiplied by renormalization factors of 0.4 to 0.7 depending on the data set, 
in particular by 0.68 in the case of our fiducial NS112 sample. We have verified that this 
procedure yields error estimates similar to those for the case when no weights are applied. 



3.3. Mean Density 

To calculate the mean galaxy density p, one can use the minimum-variance estimator 
of Davis and Huchra (1982) (also [Loveday et al. 19921 ). We can adapt it to the LCRS as 



follows. For galaxies i within redshift limits zi < zi < Z2 and absolute magnitude limits 

Ml < Mi < M2, we have 

_ _ Eiwjzi) Wi 

' Siz)wiz)i%)dz ^ ^''^ 

where S{z) is the selection function 

Siz) = / (l)oUM)dM / / 0(M)dM , (26) 

Jmax[M^i^^{z),Mi] / J Mi 



and w{z) is the weight function 



1 + fpJsSiz) 



with 



00 



J3 = 47r / r'^{r)dr {2t 







denoting the second moment of the two-point spatial correlation function ^(r). We use 
J3 ^ 10000 {h~^ Mpc)^ as found from a fit to the LCRS real space correlation function ( |Lin| 
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et al. 1996b). Recall (pobs is the Schechter function convolved with a Gaussian magnitude 
error distribution, and the weights Wi account for the sampling and completeness effects 
(2)- (4) of the previous subsection. Formally we should also incorporate the completeness 
effects (3) and (4) into the definition of S{z) above, but since these completeness corrections 
are generally not too different from unity, it is simpler to use the weights Wi instead. Also, 
for simplicity, we just use the mean sampling fraction / in the denominator of w{z) and 
ignore the field-to-field variations. 



The mean density p is calculated iteratively using equations (^) and (^). However, 
we find that for the redshift and absolute magnitude limits we are interested in, 1000 
km s^^ < cz < 60000 km s~^ and —23.0 < M < —17.5, the minimum-variance estima- 
tor gives essentially the same p as when we use inverse selection function weighting. 



w{z) = 1/S{z) 



(29) 



We will thus adopt this simpler weighting. Nevertheless, we will estimate errors 6p in the 
mean density using the minimum-variance weighting scheme 



6p 



p/f 



U:^^ Siz)wiz)i%)dz 



- 1/2 




1/2 


~ P 


.V. 





(30) 



(Here w{z) is that of equation Note also that the dependence of Sp on the actual value 
of J3 used is fairly weak.) However, this error estimate accounts just for that part due to 
galaxy clustering; we add in quadrature the uncertainty in p arising from varying M* and 
a along their joint la error ellipse. Finally, we can compute the normalization 0* in the 
Schechter function (|D by 

P 



itt,ct^'{M)dM' 

where 0' is the Schechter function with 0* set to one. As with p, the errors in 0* 
contributions both from galaxy clustering and from uncertainties in M* and a. 



(31) 
will include 



4. Luminosity Function and Mean Density of LCRS Galaxies 

Details of the STY fits to the luminosity functions of the N50, S50, N112, S112, and 
combined NS112 samples are given in Table 2a error ellipses, calculated using equation (|) 
with Ax^ = 6.17, are plotted in Figure g The STY and SWML solutions for the N112, SI 12, 
and NS112 samples are shown in Figures ^ and^. We note first that the N50 and S50 results 
seem discrepant with the NS112 result. We have verified that we can recover the correct 



M* and a parameters, using the methods of § |3l^, for the N112 and SI 12 samples culled to 
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the 50-fiber limits. It is thus unclear why the N50 and S50 samples are still discrepant; the 
trends of the discrepancies suggest some residual zero-point offset for the S50 sample and 
some remaining incompleteness in the N50 sample. Since the 50-fiber data constitute just 
20% of the whole sample, and since they are subject to larger corrections for central surface 
brightness selection (Figure |a), our final quoted results will use just the 112- fiber data, even 
though the agreement between the N112 and S112 samples is also not quite ideal. The N112 
and SI 12 2a error ellipses just touch in Figure |. The M* values do agree well, but the 
S112 has an a more positive by 0.1 than that of N112. This suggests an incompleteness at 
the faint end of the luminosity function, but examination of Figure ^ shows instead that the 
difference arises from an excess of galaxies in the S112 sample at the bright end M < —20. 
The explanation for this North/South difference is also unclear. 

We have done extensive work in eliminating residual field-to-field photometric zero-point 
offsets in our data; details of the procedures used are given in Shectman et al. (1996). In 
brief, relative zero-point offsets between our drift scans are removed by comparing data in 
overlapping regions of the scans and by comparisons to additional CCD calibration frames. 
Moreover, the effects of variable sky brightness and seeing on our isophotal magnitudes are 
dealt with by comparing, on a field-by-field basis, the offsets between galaxy isophotal and 
aperture magnitudes (using a 12"-diameter aperture that should be insensitive to seeing and 
sky brigtness changes). After the correction procedures, the mean photometric zero-points 
from slice to slice showed no offsets greater than ^ 0.015 magnitude. It is thus unlikely 
that zero-point offsets internal to the data can account for the differences seen in Figure ^ 
particularly the suggestive offset of the S50 sample. Externally, a difference in extinction 
between the North and South galactic caps comes to mind. We have not made a previous 
correction for this because 96% of our fields are at high galactic latitudes |6| > 40°. If 
we apply an extinction law Am = — 0.07(csc \b\ — 1) for |6| < 50° and Am = otherwise 
( [Sandage 1973| ), the changes in both M* and a are less than 0.02 for all four data sets. An 
additional extinction of about 0.1 magnitude in the North relative to the South could improve 
the overlap between the N112 and SI 12 error ellipses, but would not help with the 50-fiber 
results. Another possibility is that our errors are underestimated. However, we find that the 
scatter among subsamples of the four main data sets are in fact consistent with the estimated 
error ellipses for the subsamples. On the other hand, the error estimates do not account for 
any intrinsic variations in the luminosity function parameters with location in the universe. 
Thus in particular the N112/S112 discrepancy in a may hint at such cosmic variance, because 
of the large separation between the samples, but the N50/N112 and S50/S112 differences 
do not, as these last two pairs of samples are in close proximity. Another point to note is 
that the image quality for the 112-fiber photometry is appreciably improved over that of 
the early 50-fiber data photometry, due to the improved CCD's used for the 112-fiber drift 
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scans ( IShectman et al. 1995| ). In any event, whether the sample differences are due to some 
residual systematic effect(s) or not, we find that the M* values for the N50, N112, and SI 12 
samples are in good agreement, and we will see below in Table |^ that the galaxy number and 
luminosity densities for the four main data sets are also consistent within their estimated 2cj 
errors. Nevertheless, in light of the difference in the N112 and S112 a's, we will revise our 
estimate of the error in a for the NS112 sample upward from the formal value (see below). 

The likelihood-ratio and test (§ |3.1|) probabilities P listed in Table ^ show that the 
Schechter function generally gives only formally marginal fits (P > 0.001) to the luminosity 
function of the LCRS samples. Nonetheless, a visual inspection of Figures |] and ^ shows 
that the Schechter function is still a useful approximation for the absolute magnitude range 
—23.0 < M < —17.5. The SWML solutions show excesses relative to the Schechter fits at 
the faint end M > —17.5. We adopt as the LCRS luminosity function the result for the 
NS112 sample, valid for the range —23.0 < M < —17.5, with la errors: 

M* = -20.29 ±0.02 + 5 log /i 
a = -0.70 ±0.05 (32) 
(p* = 0.019 ±0.001 Mpc"^ . 

Here we conservatively estimate the uncertainty in a to cover the range in a derived from the 
N112 and S112 samples; the formal error is ±0.03. Other than this enlargement of the error 
in a, we note that the errors quoted above and elsewhere in this paper are formal values which 
should account for random uncertainties in the luminosity function parameters. Systematic 
uncertainties arising from the following sources should be kept in mind and were discussed 
previously: completeness corrections (§ p. 21) , extinction and photometric calibrations (earlier 
this section), and cosmology (§ Q). 

In our subsequent analyses of the LCRS data (apart from the present paper), we will use 
the luminosity function values in (^) for the N50, N112, and S112 data sets, but for the S50 
data set, we use its own luminosity function parameters as given in Table ||. In Figure |^ we 
plot the selection function S{cz) computed from equation (^) using the M* and a values 
for each of our four data sets and for the combined NS112 sample, while restricting the 
absolute magnitudes to —23.0 < M < —17.5. Note the typical decline at high cz caused by 
the survey's faint apparent magnitude limits. However, there is also a decline at small cz, a 
result of the survey's additional bright apparent magnitude limits. 

The mean densities computed for the LCRS samples are given in Table ^ We plot in 
Figure ^(a) the galaxy density p as a function of redshift cz, for galaxies in the NS112 sample 
within absolute magnitude limits —23.0 < M < —17.5. We also plot the cumulative density 
for galaxies lying within redshift cz, showing that the density appears to have converged 
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within the LCRS survey volume. We stop at cz = 60000 km s~^, where the selection 
function S{cz) has dropped to approximately 0.01 (Figure ^. We find a mean density p = 
0.029 ± 0.002 galaxies Mpc~^. Galaxy clustering contributes a 4% fractional uncertainty 
to p, and uncertainties in the luminosity function parameters contribute another fractional 
uncertainty of 4%. Note from Table § that despite the difference in a values, the mean 
densities of the N112 and S112 samples do in fact agree well. To compare the N112 and S112 
radial density distributions in more detail, we also plot, as a function of cz, the ratio of the 
N112 to SI 12 galaxy number densities (Figure ^jc) and galaxy number histograms (Figure [^d; 
account has been made of the differences in volume and completeness). Note that the ratios 
differ somewhat in the two figures, since the density ratios include the effect of the differing 
N112/S112 luminosity functions, but the histogram ratios do not. We find from either plot 
that for radial shells of thickness 5000 km s~^, fiuctuations in galaxy numbers and densities 
of about 20% between N112 and S112 are not uncommon. In particular, the decline seen in 
the N112 sample for cz > 40000 km s~^ is consistent with the deficit of M < M* galaxies 
observed in the N112 vs. S112 luminosity functions (Figure P). On the other hand, we also 
plot, as error bars on the ratios, the expected fiuctuations arising from the combination of 
galaxy clustering, shot noise, and uncertainties in M* and a (the last only for the density 
ratios). In computing the galaxy clustering component of the fiuctuations we have accounted 
for the geometry of the shells and used a fit to the LCRS power spectrum ( [Lin et al. 1996a ; 



this procedure is preferred over the approximate error estimate of equation |Q, which tends 
to give an overestimate given the geometry of the shells). We see that the observed density 
fiuctuations between N112 and S112 appear to be consistent with the expectations, despite 
the 2(7 discrepancies between the N112/S112 M*-a error ellipses. The N112/S112 luminosity 
functions thus do not differ significantly in normalization (as described by the galaxy number 
densities), but do differ in shape (specifically in a). 

We also tabulate in Table ^ the galaxy luminosity density pi, given by 

PL = 10-°-^^'XM)rfM . (33) 

J-23.0 

(Letting the integration limits go to infinity only increases pi by 3%.) The NS112 value 
of (1.4 =b 0.1) X 10^ Lq h Mpc~^, when increased by 25% as a rough isophotal-to-total 
light correction, implies a critical mass-to-light ratio {M/L) ^ 1.6 x 10'^ {M/L)q to close 
the universe, in agreement with previous determinations (e.g., [Kirshner et al. 1983|; EEP; 



Marzke et al. 1994a| ). Alternatively, we obtain a cosmological density parameter Vt 



0.2 [(M/L) / 300 h {M/L)q], if we set M/L = 300 h {M/L)q, a typical virial mass- 
to-light ratio estimated from clusters and groups (e.g., a recent determination from the the 
CNOC cluster sample gives M/L = 295 h {M/L)q in Gunn r; [Carlberg et al. 1996| ). 

In addition, we have also derived 0* for the NS112 sample by normalizing to our galaxy 
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number count distribution dn{m). We use 

m+Am roo rM+AM 

dn{m) = / dz{dV/dz) / ^obs{M')dM' , (34) 

m Jo Jm 

and do a x^-fit to find 0*, assuming Poisson errors on the counts, and fixing M* and a at their 
maximum hkehhood values. For the counts we consider all objects in our CCD photometric 
catalog with isophotal magnitudes 15.0 < m < 17.8, except that we make a 3% correction 
for stellar contamination, as inferred from the fraction of stars in our spectroscopic sample. 
We find (p* = 0.021 Mpc~^, but to properly compare against our redshift survey value, 
we need to make a 4% reduction to account for the low surface brightness galaxies excluded 
from the redshift survey. We get 0* = 0.020 Mpc~'^, in good agreement with the value 
(p* = 0.019 ± 0.001 Mpc~^ derived from the redshift survey. We show in Figure ^ the 
dn{m) distribution for the NS112 sample and the prediction based on 0* = 0.021 Mpc~^, 
M* = -20.29, and a = -0.70. 

We have also separated the NS112 sample into two groups by [Oil] 3727 emission, 
defining emission galaxies to be those with 3727 equivalent widths Wx > 5 A. The emission 
galaxies constitute 40% of the total sample. The fit details are given in Table 0, the error 
ellipses are plotted in Figure |TU], and the SWML and STY solutions are shown in Figure |TT[ 



Note that for the S112 data, we have used separate surface brightness completeness functions 
G{M), shown in Figure 0(b), for the emission and non-emission galaxies. We see that 
emission and non-emission galaxies have very different luminosity functions, with emission 
galaxies dominating the faint end and non-emission galaxies prevailing in the bright end. 
In particular, emission galaxies constitute about 15% of the total at M = —22, rising to 
about 80% at M = —17. The two populations show about equal number densities near 
the full-sample M*. The emission galaxies show a much steeper faint-end slope a (—0.9 vs. 
—0.3) and a slightly fainter M* (—20.0 vs. —20.2) relative to the non-emission galaxies. Also, 
note from Figure |TT| that the excesses at M > —17.5 relative to the Schechter fits are less 
pronounced than they were for the full sample, though for consistency we use the same fit 
range —23.0 < M < —17.5. Moreover, the probabilities in Table |^ indicate that a Schechter 
function is a good description of the emission galaxy luminosity function, though again visual 
inspection of Figure indicates that the Schechter function provides reasonable (though 
not always formally good) approximations to the luminosity functions of both subsamples. 
In figure |^(b), we plot the space density of emission and non-emission galaxies. Finally, we 
note that instead of our formal la errors for a listed in Table ^ we quote instead an error 
of ±0.1 for both the emission and non-emission galaxies. We do this because of our earlier 
N112/S112 difference in a, and because of another effect to be discussed at the end of § 0. 
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5. Comparisons to Other Surveys 

We next compare our luminosity function to that from three other large, local, and 
optically-selected redshift surveys. The luminosity functions from the Stromlo-APM (jLove- 



|day et al. 1992 ), Center for Astrophysics (CfA; [Marzke et al. 1994a ), and second Southern 



Sky (SSRS2; |da Costa et al. 1994|) redshift surveys are detailed in Table || and are plotted 



against that of the LCRS NS112 sample in Figure |T^. Because these other surveys were 
selected in the blue, while the LCRS was selected in the red, for the comparison we need to 
shift the zero-points of the various magnitude systems. (Shifting the zero-points also roughly 
accounts for other differences in the magnitude systems; e.g., the fact that LCRS magnitudes 
are explicitly isophotal whereas APM magnitudes have an isophotal-to-total magnitude cor- 
rection applied. However, applying simple constant offsets does neglect more complicated 
effects, due for example to galaxy colors and morphologies. Be aware that the specific dif- 
ferences we find below in the luminosity functions will necessarily be sensitive to the exact 
shifts in absolute magnitudes we apply.) For Stromlo-APM we have offset its Mi,j magnitude 
scale by 1.1 magnitude to approximately match the mean rest-frame color {hj — R)o = 1-1 
of LCRS galaxies ( rTucker 199^ P^ucker et al.^995| ). We shift the SSRS2 Mb{o) magnitude 



scale by the same amount, since the SSRS2 and bj zero-points agree within 0.2 magnitude 
( |da Costa et al. 1994|) . However, we shift the Zwicky Mz magnitude scale for the CfA survey 



by another 0.7 magnitude to match the CfA and Stromlo-APM M* values. Note first from 
the figure that the Stromlo-APM and LCRS results show consistent normalizations and are 
fairly similar in overall shape, despite differences in their respective Schechter parameters, 
particularly in the faint-end slope a. There are instead differences, for Mr < —21.5 and 
Mj. > —19, which are responsible for the shallow faint-end slope a = —0.7 found in the 
LCRS, compared to the more typical flat value a ~ — 1 seen in the Stromlo-APM and other 
previous surveys. We emphasize that even though the LCRS is slightly biased against in- 
trinsically faint galaxies (because of the survey selection criteria and the correlation between 
luminosity and central surface brightness), the gross features of our luminosity function are 
quite similar to an a = —1 Schechter function. However, such a Schechter function does 
not provide as good a fit in detail to the LCRS data. As another aid in the comparison, 
note that for -1 < a < -0.7, we roughly have (Apl/pl) ~ (A070*) - 0.9AM* + |(Aa/a) 
(good to about 0.1 in Ap^^/Pi)- Thus the LCRS and Stromlo-APM surveys have about the 
same luminosity density pi despite their different Schechter parameters (Table 1^): the 30% 
difference in 0* is more than compensated by the 30% difference in M* (after noting the 1.1 
magnitude zero-point shift) and one third of the 30% difference in a. 

We note next from Figure |12| that the CfA and SSRS2 luminosity functions show some 
conspicuous differences compared to the LCRS and Stromlo-APM results. The CfA luminos- 
ity function clearly shows a higher normalization, though its shape is similar as it has a = —1. 
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on the other hand, the SSRS2 luminosity function is quite similar to those of the LCRS and 
Stromlo-APM surveys for M^j < —19, but for fainter magnitudes its a = —1.2 makes the 
SSRS2 result overshoot the LCRS and Stromlo-APM luminosity functions. Though not 
shown in the figure, the CfA luminosity function also has a similar excess at faint magni- 



tudes Mz > — 16 relative to an extrapolation of the a = —1 fit (see Figure 9 of [Marzke 



et al. 1994a|) . It is unclear why such a steep faint end (see below) is not observed in the 
LCRS or Stromlo-APM surveys. However, the high normalization of the CfA survey may 
be explained as the result of local density inhomogeneities (e.g., large coherent structures 
like the "Great Wall"; [Geller fc Huchra 1989|) present in the CfA survey volume. Both the 
CfA and SSRS2 surveys probe a smaller and more local volume of space compared to either 
the LCRS or Stromlo-APM surveys; the mean survey depths are about 7500 km s~^ for 
both CfA and SSRS2, 15000 km s'^ for Stromlo-APM, and 30000 km s"^ for the LCRS. In 
fact Marzke et al. (1994a) and da Costa et al. (1994) both point out that their respective 
surveys are not yet large enough to be "fair" samples of the universe. Both the LCRS and 
Stromlo-APM surveys should come closer to representing fair samples, and it is remarkable 
that their respective luminosity functions are quite similar in both normalization and shape. 

We do not have morphological classifications for our survey galaxies. Nonetheless, 
though there is not a one-to-one correspondence between 3727 emission and morphological 
type, we do expect our emission galaxies to tend to be late-type spirals and our non-emission 
galaxies to be early-type ellipticals and SO's. We also tabulate in Table ^ the luminosity func- 
tion for galaxies of different morphological types computed for the Stromlo-APM and CfA 
( [Marzke et al. 1994b| ) surveys. The sharp distinction in a between the early- and late-type 



luminosity functions in the Stromlo-APM survey {a = +0.2 and —0.8 respectively) quali- 
tatively agrees with the emission vs. non-emission differences seen in the LCRS. However, 
Zucca et al. (1994) and Marzke et al. (1994b) suggest that incomplete classification of early- 
type galaxies may be responsible for the difference seen in the Stromlo-APM survey. In 
fact, there is no significant difference in the luminosity functions of E-SO and spiral (Sa-Sd) 
luminosity functions in the CfA survey {a = —0.8 to —0.9). However, Marzke et al. (1994b) 
do find that their faint end is dominated by Magellanic spiral and irregular galaxies, with 
a very steep a = —1.87. Moreover, da Costa et al. (1994) attribute their steep a = —1.2 
to the dominance of blue galaxies at the faint end. These CfA and SSRS2 results are in 
qualitative agreement with the LCRS findings, as we expect such late-type blue galaxies to 
be gas-rich and to exhibit strong emission lines. 
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6. Relation to Galaxy Counts and Galaxy Evolution 



The problem of the excess faint blue galaxy counts over the expectations of the simplest 

(see reviews by [Koo fc| 



no-evolution models has continued to be an area of intense research 



Kron 1992 ; Koo 1995 ). The absence in deep pencil beam redshift surveys ([Broadhurst et 



al. 1988| ; Polless et al. 1990| , 1993; plazebrook et al. 1995a|) of high redshift z > 1 galaxies, 
which are predicted by simple luminosity evolution models for the excess counts, has led 
to proposals of more exotic explanations. For example, the excess blue galaxies may be 
star-bursting dwarfs which have since faded by the present day (e.g., [Broadhurst et al. 1988 



Cowie et al. 1991] [Babul fc Rees 19"9^ ), or alternatively, extensive galaxy merging may have 
occured to remove the excess galaxy population by 2; = ([Broadhurst et al. 1992[ ; [Rocca- 
Volmerange fc Guiderdoni 1990[ ). 



In contrast, other workers ([Gronwall fc Koo 1995[ ; [Koo et al. 1993[ ; [Metcalfe et al. 1991 



1995b) have argued that the data may still be reasonably fit by mild evolution models, given 
the uncertainties in the faint end slope and in the normalization of the local luminosity 
function; note in particular the steep a of the Sm-Im sample of Marzke et al. (1994b), and 
the difference in normalization between the Stromlo-APM and CfA surveys. In fact, even 
before the CfA results, the counts in the APM galaxy catalog ( [Maddox et al. 199Qa[ ) at bright 
magnitudes bj = 20.5 showed a factor of 2 excess relative to no-evolution predictions based 
on the counts at bj ^ 17 ( [Maddox et al. 1990 b[ ). Either significant evolution occurred at low 
redshifts z ^ 0.1 ( [Maddox et al. 1990b[) , or there are large-scale structures (voids) and/or 
magnitude scale errors ( [Metcalfe et al. 1995a[ ) in the APM survey which are biasing the 
6j ~ 17 counts systematically low. Choosing to normalize to the galaxy counts at bj ~ 20 
corresponds to using 0* = 0.03 Mpc~^, twice the Stromlo-APM value, and using this high 
normalization reduces the amount of evolution needed to reconcile with the faint count data. 

There is now evidence from deep redshift surveys for evolution in the luminosity function 
since redshifts z ~ 0.5 ( [ti^ales 1993| ; Lonsdale & Uhokshi 1993), but one which depends on 
galaxy type. Blue, star-forming galaxies appear to exhibit an increased density and/or 
steepening faint-end slope with increasing redshift, whereas the red galaxy population shows 
no evolution with z ([Lilly et al. 1995] ; Polless 1995[) . There is also evidence of an increase 
in star formation rate with redshift and the association of merging activity with increased 
star formation ( PoUess 1995[ ; [Colless et al. 199^ ). In particular, the number of star- forming 
galaxies, defined as those with [Oil] 3727 W\ > 20 A, appears to be 5-10 times greater at 
z > 0.5 than at z = 0, for absolute magnitudes below and including M* ( [Colless 1995[ ). 
In addition, recent galaxy number counts from the Hubble Space Telescope (HST) Medium 
Deep Survey ( [Glazebrook et al. 1995 b[ ; Driver et al. 1995a,b) indicate that at -B ~ 24, 
the counts are dominated by galaxies with irregular and peculiar morphologies, which show 
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increased numbers relative to their proportions in local surveys. However, the elliptical and 
spiral counts are as expected with no evolution needed, if one adopts the Stromlo-APM 
luminosity function, but with the high normalization 0* = 0.03 Mpc~^. 

It is beyond the scope of this paper to undertake detailed modelling of the evolution of 
galaxies, but results from the local LCRS sample bear on two points which may be connected 
to understanding the faint blue galaxies. One is the normalization of the local luminosity 
function and the other is the rate of evolution for galaxies with active star formation. First, 
we note that our luminosity function normalization is consistent with the low value from 
the Stromlo-APM survey, not with the high value implied by the bj ~ 20 APM galaxy 
counts. Thus results arguing for no or mild evolution, e.g. the luminosity function models 
of Gronwall & Koo (1995) or the HST elliptical and spiral counts, will need to be modified 
if our low normalization is not an underestimate. Because of our large survey depth, up 
to z = 0.2, large-scale structures should not bias our results low (cf. Figure ^). The 
normalizations for the two completely independent samples N112 and S112, which probe 
opposite directions in the sky, agree well with each other. Also, the remaining 4% of low 
central surface brightness galaxies that are not part of the redshift survey will not lead to 
such a large underestimate. We saw before that our 0* normalization is consistent with our 
total number counts, including the low surface brightness galaxies. Any potential factor of 
2 underestimate will have to come from galaxies completely missing from our original CCD 
photometric catalog, because they lie below our effective surface brightness limit. Such a 
bias against very low surface brightness galaxies, as has been argued for by, e.g., McGaugh 
(1994) and Ferguson & McGaugh (1995), would be a potential problem for all bright galaxy 
surveys, not just the LCRS. We have derived a firm value for the density normalization of 
galaxies with central surface brightnesses above our explicit cut line(s), examples of which 
are shown in Figure |I]. A final observation here is that our galaxy counts, for our hybrid 
Kron-Cousins i? = 15 — 17.8, correspond roughly to 6j = 16 — 19. However, whereas the 
APM counts over this range seem too steep at face value to be consistent with no evolution, 
even at such bright magnitudes ( [Maddox et al. 1990 b| ), the Las Campanas counts shown in 



Figure |^ at comparable magnitudes suggest no such evidence for evolution at low redshifts 
z ^ 0.1. 

Second, we address the possibility of an increase in the number of star-forming galxies at 
redshifts z < 0.2 explored by the LCRS, as has been observed by redshifts z ^ 0.3 in deeper 
surveys ( [Colless 1995| ), by checking whether the average equivalent width of [OH] 3727 varies 



with redshift. To factor out the dependence of [OH] 3727 Wx on galaxy luminosity, we plot 



in Figure |T3| the average Wx as a function of z for seven separate absolute magnitude bins, 
each one magnitude wide, from M = —16 to —23. We use only our NS112 data. Our results 
indicate that the average equivalent widths generally do not change with redshift, although 
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for galaxies fainter than M = —19, the redshift range is not large. A similar conclusion 
may also be drawn by plotting the fraction of galaxies with, say, W\ > 5 A or 20 A. In 
light of the error bars, the most significant trend seen is that for the —20 > M > —21 
bin, specifically a rise of about 1 A in the average Wx at z ~ 0.1. This may be due to 
large-scale structure, as the same trend is less apparent in the N112 sample than in the S112 
sample. Alternatively, it may be an aperture effect, in that our 3.5"-diameter fibers cover 
5 h^^ kpc on a galaxy aX z = 0.1, but only half as much at z = 0.05, so that there would 
be a bias toward lower equivalent widths at lower redshifts if less emission arose from the 
centers of galaxies compared to that from the outer regions. An argument against the latter 
explanation is that we do not see any consistent systematic trend of this kind in the other 
absolute magnitude bins. In any event, the evidence for changes in Wx is not dramatic in 
any of the luminosity bins. Since 3727 equivalent width is a surrogate for the star-formation 
rate ( [Kennicutt 1992| ), we conclude that there is no evidence for significant evolution in the 
star-forming properties of LCRS galaxies over the redshift range we sample. We emphasize 
that this conclusion is only valid over the respective redshift ranges observable for galaxies 
of different luminosities, and that galaxies fainter than M* + 1 ^ —19 are not seen above 
z = 0.1. For galaxies brighter than M*, though, our no-evolution result does extend up to 
z = 0.2. However, note that our null results say nothing about the above-mentioned evidence 
of increasing numbers of star-forming galaxies at higher redshifts than that constrained by 
the LCRS, specifically an increase for sub-M* galaxies by 2; ~ 0.3 and for M* galaxies by 
z ^ 0.5 (ICoUess 19951 ). 

A final note is that the small increase in Wx seen for M* galaxies does have some 
impact on the a value derived earlier for our emission and non-emission galaxy samples. In 
fact, in our earlier derivation we actually applied redshift-dependent weights, analogous to 
the completeness corrections discussed in § |3.2| , to account for any changes in the emission 
(or non-emission) galaxy fraction as a function of z. The weights are calculated separately 
for galaxies of different absolute luminosities. If we had not applied such corrections, the 
maximum-likelihood solution for the emission galaxy a would steepen by —0.1, while the non- 
emission a would flatten by -1-0.1. Leaving out the additional corrections leads primarily to 
a systematic increase in the emission galaxy density with redshift (because of the steeper a), 
so that the sum of the cumulative emission and non-emission galaxy densities overshoots that 
for all galaxies by about 25% at the cz = 60000 km s~^ limit. In light of the lack of evolution 
seen in the Figure 13 results, which do not depend on the luminosity function, we feel justified 
in applying the z-dependent corrections, resulting in more reasonable agreement between the 
total galaxy density and the sum of the emission and non-emission galaxy densities given in 
Table 0. We increase our estimated errors in the emission and non-emission a's to ±0.1 in 
light of the above effects, as as we alluded to before at the end of § ||. 
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7. Conclusions 

We have used inhomogeneity-independent maximum-likelihood methods to calculate the 
luminosity function of LCRS galaxies in the (hybrid) Kron-Cousins R absolute magnitude 
range —23.0 < M — 51og/i < —17.5, where we find that it may be reasonably described 
by a Schechter function with M* = -20.29 ± 0.02 + 51og/i, a = -0.70 ± 0.05, and 0* = 
0.019 ± 0.001 Mpc~^. There appears to be an excess relative to the Schechter function 
at the faint end M > —17.5. The mean galaxy density for the range —23.0 < M < —17.5 is 
0.029 ±0.002 galaxies Mpc"^ out to a redshift of 60000 km s'^, and the run of cumulative 
density with redshift suggests that the density has converged within the LCRS survey volume. 
The measured (isophotal) galaxy luminosity density is = (1.4 ± 0.1) x 10^ Lq h Mpc^^, 
resulting in ^ 0.2, for a typical virial mass-to-light ratio of 300 h {M/L)q and a rough 
isophotal-to-total light correction. We have also divided our sample by the cut [Oil] 3727 
equivalent width Wx > 5 A , and showed that the emission and non-emission galaxies so 
defined have significantly different luminosity functions, with emission galaxies dominating 
the faint end {a — —0.9 ± 0.1), non-emission galaxies dominating the bright end {a — 
—0.3 ± 0.1), and about equal numbers near the full-sample M*. 

We then compared our results to those from the optically-selected Stromlo-APM, CfA, 
and SSRS2 redshift surveys. Despite the fact that the Stromlo-APM survey has an a f=:i — 1, 
and the LCRS has a — —0.7, the shapes of two luminosity functions are remarkably similar, 
though they do differ in specific details. Moreover, our luminosity function normalization 
is consistent with that of Stromlo-APM; the CfA normalization is higher but this may be 
biased by density inhomogeneities in its smaller and more local survey volume. Finally, our 
emission galaxy results confirm the picture seen in these previous surveys that the faint end 
of the luminosity function is dominated by late-type, blue galaxies which are presumably 
also gas-rich and strong in [Oil] 3727 emission. 

Finally, two implications of our results may be made with respect to the faint galaxy 
counts and to galaxy evolution. First, our normalization of the luminosity function is a 
factor of two below the "high" normalization often adopted from the hj ~ 20 APM galaxy 
counts in calculations used to match fainter galaxy count data. If our normalization is not 
an underestimate (e.g., due to surface brightness effects), then it implies that more evolution 
is required in the models, compared to using the high APM count normalization. Second, for 
the redshifts sampled in the Las Campanas survey, up to about z — 0.2 for galaxies brighter 
than M*, our results do not indicate a significant change of average [Oil] 3727 equivalent 
width with redshift. The star-formation rate, as measured by [Oil] 3727, has apparently not 
evolved rapidly over the volume sampled by the LCRS, unlike the large increases seen in the 
numbers of blue star-forming galaxies in smaller, deeper surveys at ^ > 0.3. 
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In this paper we have used the Las Campanas Redshift Survey to provide an accurate 
determination of the local luminosity function for the largest galaxy sample thus analyzed. 
Our red i?-band results are also complementary to previous results in blue photometric 
bands. In addition, we have begun to explore the luminosity function for galaxy samples 
divided by emission properties, and in future analyses we will examine galaxy samples defined 
by &j — -R colors (which are available for our Southern galaxies from APM observations), 
and if possible, also by basic morphological information to be determined from our drift- 
scan images. The Las Campanas Redshift Survey should thereby provide a useful and 
large database of local galaxy properties, including luminosity, emission, color, and possibly 
morphological information, upon which to anchor models of faint galaxy counts and galaxy 
evolution at higher redshifts. 

We thank Ron Marzke for useful discussions and the anonymous referee for helpful 
suggestions. The Las Campanas Redshift Survey has been supported by NSF grants AST 
87-17207, AST 89-21326, and AST 92-20460. HL also acknowledges support from NASA 
grant NGT-51093. 

Note added in proof: the difference noted here between the luminosity functions of 
emission and nonemission galaxies is similar to that found in preliminary analysis of the 
ESP redshift survey galaxies, e.g., Vettolani et al. (1994), lAU Symposium No. 161, 687. 
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Table 1. LCRS Subsamples 



Sample 


N " 




mi 


1712 


% LSB d 


North 50-fiber (N50) 


2706 


0.66 


16.0 


17.3 


22% 


South 50-fiber (S50) 


2054 


0.50 


16.0 


17.3 


21% 


North 112-fiber (N112) 


8552 


0.71 


15.0 


17.7 


4% 


South 112-fiber (S112) 


10378 


0.69 


15.0 


17.7 


9% 



^Number of galaxies observed which lie within the survey photometric limits and 
geometric borders. 

'^Mean sampling fraction for galaxies within photometric limits and geometric 
borders. 

■^Nominal isophotal magnitude limits: rui <m < m^.- 

'^Percentage of objects excluded for low central surface brightness but otherwise 
within the isophotal limits. 
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Table 2. Luminosity Function Parameters 



Sample 


N 


M*-51og/i^ 


a 








LCRS 














N50 


2690 


-20.24 ±0.07 


-0.45 ±0.11 


0.017 ±0.003 


0.20 


0.71 


S50 


2043 


-20.55 ±0.10 


-0.74 ±0.13 


0.016 ±0.003 


0.03 


0.04 


N112 


8399 


-20.28 ±0.03 


-0.75 ±0.04 


0.018 ±0.002 


0.002 


0.004 






on oo 4- n DQ 


— U.OO ± U.U4 


U.UZU ± U.UUZ 


U.Ui 


U.UO 


NS112 


18678 


-20.29 ±0.02 


-0.70 ±0.03 


0.019 ±0.001 


0.001 


0.009 


3727 > 5 A 


7312 


-20.03 ±0.03 


-0.90 ±0.04 


0.013 ±0.001 


0.24 


0.45 


3727 Wx<5A 


11366 


-20.22 ±0.02 


-0.27 ±0.04 


0.011 ±0.001 


0.002 


0.02 


Dtromio-Ar iVi 














All 


1658 


-19.50 ±0.13 


-0.97 ±0.15 


0.0140 ±0.0017 


0.65 




Early 


311 


-19.71 ± 0.25 


±0.20 ± 0.35 




0.64 




Late 


999 


-19.40 ±0.16 


-0.80 ±0.20 




0.87 




CfA ^ 














All 


9063 


-18.8 ±0.3 


-1.0 ±0.2 


0.04 ±0.01 


0.25 




E-SO 




-18.87 


-0.92 


0.010 






Sa-Sd 




-18.76 


-0.81 


0.015 






Sm-Im 




-18.79 


-1.87 


0.0006 ± 0.0002 


0.46 




SSRS2 ^ 


2919 


-19.50 


-1.20 


0.015 ±0.003 







''Ho = 100 h km s"^ Mpc'^ 

^Pi is probability for likelihood ratio test; P2 is probability for dN{M) test. See § 3.1. 

'^Magnitudes are in our "hybrid" Kron-Cousins R system and absolute magnitudes are 
restricted to —23.0 < M — 51og/i < —17.5. 0* is calculated for redshifts 1000 km s~^ < 
cz < 60000 km s~^. Errors are formal la uncertainties; but also see discussion in §§ 4 and 
6. 

■^Loveday et al. 1992. Magnitudes in bj system. 
^Marzke et al. 1994a,b. Magnitudes in Zwicky system, 
^da Costa et al. 1994. Magnitudes in -B(O) system. 
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Table 3. Number and Luminosity Densities 



Sample 




PL(lO^L0/iMpc-=^) ^'^ 


N50 


0.020 ±0.004 


1.2 ±0.1 


S50 


0.026 ±0.006 


1.5 ±0.2 


N112 


0.028 ±0.002 


1.3 ±0.1 


S112 


0.029 ±0.002 


1.5 ±0.1 


NS112 


0.029 ±0.002 


1.4 ±0.1 


3727 Wx>5A 


0.021 ±0.002 


0.74 ± 0.05 


3727 W^A < 5 A 


0.011 ±0.001 


0.78 ± 0.04 



^We take Hq = 100 h km s"^ Mpc~^. p and pi arc calculated for 
absolute magnitudes —23.0 < M — Slog/?. < —17.5 and redshifts 1000 
km s~^ < cz < 60000 km s~^. Errors are formal la uncertainties. 

*^We assume Rq = 4.52 (M. Pinsonneault, private communication). 
Also, be aware that our magnitudes are isophotal and are not total; we 
estimate that our isophotal magnitudes miss roughly 25% of the total 
light on average. 
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-12° Slice, Epoch 1992.3 



"1 I I ^ I I I I ^ I I I I ^ I I I I 
example 1 12-fiber limits 

1 5.0 < m < 1 7.7, m^ < m^ = 1 9.35 - 0.5(1 7.7-m) 
example 50-fiber limits 

16.0 < m < 17.3, m^ < m^.^,, = 18.65 - 0.5(1 7.3-m) 
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16 16.5 17 

isopliotal magnitude m 
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Fig. 1. — Examples of the photometric selection criteria in the isophotal magnitude (m) - 
central magnitude (mc) plane apphed to the 112- and 50-fiber LCRS data. 




Fig. 2. — (a-b) The identification rate (galaxies plus stars) Fi{m) as a function of apparent 
magnitude m for the four main subsets of LCRS data. Vertical lines show the nominal 
isophotal magnitude limits, (c-d) The identification rate ^2(777.0 — mc^nmit) as a function of 
"distance" rric — mc^umit from the central magnitude limit. 
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Fig. 3. — (a) Completeness function G{M) vs. absolute magnitude M. G{M) is calculated 
from the North 112 sample by culling it to match the central surface brightness selection 
criteria of the North 50/South 50 and South 112 samples. Horizontal lines show the expected 
values of G{M) if the completeness does not vary with absolute magnitude, (b) The North 
112 sample is again culled to the South 112 sample selection criteria, but now G{M) is 
separately plotted for all galaxies and for subsets divided by [Oil] 3727 equivalent width. 
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Fig. 4. — 2a error ellipses in M* and a for the STY maximum likelihood fits to the four 
main subsets of the LCRS data and to the combined North+South 112 sample. Details of 
the fits are given in Table 2. 
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absolute magnitude M 



Fig. 5. — Luminosity function for the combined North+South 112 sample. The squares show 
the SWML solution, plotted with la errors, and the line shows the STY solution. Details 
of the STY fit is given in Table 2. The dotted line shows the extrapolation of the fit to 
M > -17.5. 
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absolute magnitude M 



Fig. 6. — Luminosity functions for the North 112 and South 112 samples. The points show 
the SWML sohitions, plotted with la errors. The lines show the STY solutions, with details 
of the fits given in Table 2. Dotted lines show the extrapolation of the fits to M > —17.5. 
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Fig. 7. — The selection function S{cz) , computed using the luminosity function parameters 
and nominal apparent magnitude limits of each of the North 50, South 50, North 112, South 
112, and combined North+South 112 samples. Absolute magnitude limits —23.0 < M < 
—17.5 have also been applied. 
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Fig. 8. — The galaxy density p as a function of redshift cz for: (a) the whole North+South 
112 sample; and (b) the North+South 112 sample divided by 3727 emission. The horizontal 
hues show the respective mean densities; see Table 3. The ratio of the North 112 over the 
South 112 (c) galaxy number densities and (d) normalized galaxy number histograms are 
also shown. The errors bars in (c-d) are described in the text. In all cases, the absolute 
magnitude range is —23.0 < M < —17.5. 
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- I 



15 15.5 16 16.5 17 17.5 

isophotal magnitude m 



Fig. 9. — The number count n as a function of isophotal magnitude m for the North+South 
112 photometric catalog. Poisson errors are shown. A 3% correction has been apphed to 
account for stellar contamination. The solid line shows the prediction derived from the 
NS112 luminosity function, with normalization fixed by the number counts. 
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Fig. 10. — 2a error ellipses in M* and a for the STY maximum likelihood fits to the combined 
North+South 112 sample, whole and divided into 2 subsets by [Oil] 3727 equivalent width. 
Details of the fits are given in Table 2. 
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Fig. 11. — Luminosity functions for the combined North+South 112 sample, divided by [Oil] 
3727 equivalent width. The points show the SWML solutions, plotted with la errors. The 
hues show the STY solutions, with details of the fits given in Table 2. Dotted lines show the 
extrapolation of the fits to M > —17.5. 
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blue absolute magnitudes Mj,j, M^.p,, M,-0.7 




-23 -22 -21 -20 -19 -18 -17 -16 

red absolute magnitude 



Fig. 12. — The luminosity function for the combined North+South 112 LCRS sample, com- 
pared to those of the blue-selected Stromlo-APM, CfA, and SSRS2 surveys. Note that the 
blue Mbj and Mb{o) (top) and red Mr (bottom) absolute magnitude scales are offset in or- 
der to match the average color {bj — R)o = 1.1 of LCRS galaxies. The blue Zwicky Mz 
magnitude scale has been shifted to match the CfA and Stromlo-APM M* values. The CfA 
results show an excess relative to the plotted fit at Mz > —16; see Marzke et al. (1994a). 
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Fig. 13. — The average [Oil] 3727 equivalent width as a function of redshift, computed for 
galaxies in the LCRS North+South 112 sample. The galaxies have been divided into seven 
bins in absolute magnitude, each one magnitude wide, ranging from M = — 16 to —23. The 
errors shown are standard deviations of the mean. 



